Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 3 March 2013 (MN style file v2.2) 



The Lives of High Redshift Mergers 



O 



(N 

O: 

u 

o 



> 
m 

o 

(N 



Tom McCavana^*, Miroslav Micic^'^, Geraint F. Lewis^, Manodeep Sinha^, 
Sanjib Sharma^, Kelly Holley-Bockelmann^ and Joss Bland-Hawthorn^ 

^Sydney Institute for Astronomy, School of Physics, A28, The University of Sydney, NSW 2006, Australia 
^Astronomical Observatory Belgrade, Volgina 7, P.O.Box 74 11060 Belgrade, Serbia 
'^Department of Physics & Astronomy, Vanderbilt University, Nashville, TN, USA 



3 March 2013 



ABSTRACT 

We present a comparative study of recent works on merger-timescales with dy- 
namical friction and find a strong contrast between idealized / isolat ed mergers 
(iBovlan-Kolchin et al.ir2008| ) and mergers from a cosmological volume (jJiang et al.l 
|2008| ). Our study measures the duration of mergers in a cosmological N-body simula- 
tion of dark matter, with emphasis on higher redshifts (z^lO) and a lower mass range. 
In our analysis we consider and compare two merger definitions; t idal disrup tion and 
coalescence. We find that the merger-time formula proposed by iJiang etal.l (j2068l) 
describes our results well and conclude that cosmologically motivated merger-time 
formulae provide a more versatile and statistically robust approximation for practical 
applications such as semi-analytic/hybrid models. 

Key words: (cosmology:) dark matter, methods: n-body simulations, galaxies: evo- 
lution, galaxies: kinematics and dynamics 



1 INTRODUCTION 

Galaxy mergers are aided by dynamical friction, a process 
that extracts energy and angular momentum from the orbit 
of an incoming galaxy through a gravitation al wake that acts 
as a drag force. Its theoretical conception l|Chandrasekharl 
1 19431 ) was based on a point mass traveling through an 
in-finite, uniform and non-self gravitating medium. With 
relatively few modifications to account for mass loss in 
the satellite galaxy, it is used to estimate the duration of 
galaxy mergers. As such it is a centrepiece of galaxy evo- 
lution; for example, in semi-analytical and hybrid models 
of galaxy evolution the accuracy of analytical merger-times 
is directly responsible for quantities such as stellar mass, 
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ies, known as proto-galaxies, are born at the heart of these 
DMHs from gas that has been virialized, radiatively cooled 
and collapsed into a centrifugally supported disk. These 
proto-galaxies are the birth place of quasars and AGNs 
which later evolve into the spiral and elliptical galaxies ob- 
served today. As dark matter halos merge into larger struc- 
tures so do the central proto-galaxies. This process (thought 
to be responsible for the AGN duty cycle) feeds the central 
AGN by providing a reservoir of gas which is then accreted 
onto the central massive black hole. 

The dynamics of DMH mergers is described with rea- 
sonable success by a combination of gravitational free-fall, 
dynamical friction and tidal stripping. Chandrasekhars dy- 
namical friction has been analytically developed to predict 
the coa lescence time for a circul ar orbit from some initial 
radius (|Binnev fc Tremainelll987l ). Most versions of the ap- 
proximating formula follow the form: 



In current ACDM cosmology the first objects to form 
in the universe we re made out of dark matter dZwicky 193^; 
iRubin et aPllQSGl ). Widely known as Dark Matter Halos 
(DMHs hereafter), their formation through violent collapse 
and relaxation results in a shallow inner density profile, with 
the remainder of accreted d ark matter trailing behind in a 
steep outer slope (jCiardi fc Ferrarai2005i ). The first galax- 
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= Cdf 



In A Vc{rvir) ' 



(1) 



where Cdf is a constant, Mh and Ma are the masses of the 
host and satellite DMHs respectfully a nd ,f (e) is a fun ction of 
circularity introduced analytically bv lLacev fc Colj f|1993,). 
InA is the Coulomb logarithm (discussed below), Vc{rvir) 
is the circular velocity at the virial radius (r„ir) of the host 
halo, and r« is a radius that varies depending on assumptions 
made. For bound circular orbits, r* is the initial radius 
(,Binnev fc Tremaind [l987l ). In the more general case of a 
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boun d orbit (radial or tangential) r* is rc{E) (|Lacev fc Cold 
is the radius that a circular orbit would have 
with the same orbital energy, r, is also often approximated 
by r^ir as most works assume r; ~ r^ir or rc{E)/rvir ~ 1- 

Studies to improve the approximation of DMH 
merger-times (and by proxy the effects of dynamical 
friction) are numerous and range from analytical/ semi- 



analytical (e.glWeinberg 19891: Tremaine fc Weinberg 1984 
Lacev fc Col? Il993l: iBenson et all l2005l: iTavlor fc Babull 



200ll : Ivan den Bo sch et alj Il999ir jcolpi et alj 1 19991) 



numerical (e.g. Ahmad &: Cohen _ 1973 : Lin fc Tremaine 



to 



19831 : iTormen et aLlll998l : iHashimoto et al.H2003l : iFuiii et al 



20061 ) ■ Significant work has also been done to assess the ac- 



curacy of the Coulomb logarithm (see e.g. Jiang et al. 2008) . 
Analytically, A is found to be a function of the maximum 
and minimum impact parameters. For extended bodies a 
better agreement is found with InA — ln(l -I- Mh/Ms), where 
the satellite mass is determined at in-fall. 

In order to refine the accuracy of dynamical friction 
formulae, Boylan-Kolchin et al. (2008, BKQ8 hereafter) ex- 
plore merger timescales over a range of mass ratios, circu- 
larities and orbital energies using a suite of simulations of 
isolated mergers. In their analysis they build two well re- 
solved Hernquist halos and merge them in various combi- 
nations of initial orbital parameters. They also explore the 
infiuence of baryons and find that including a stellar bulge 
{Mt,/MDM = 0.05) shortens the merger-time by ^ 10%. 
BK08 produce a empirical formula that relates merger-time 
to a mergers initial orbital parameters. They argue that a 
typical application of equation [T] systematically underesti- 
mates the merger-times seen in their simulation. This dis- 
crepancy grows as the host becomes larger than the satellite. 

As one of the first simulations to study dynamical fric- 
tion in a cosmological context, Jiang et al. (2008, J08 here- 
after) developed a merger-time formula using direct mea- 
surements from a high resolution cosmological SPH /N-body 
simulation. Their measurements provide a more physically 
motivated approach for estimating merger-times. Galaxies 
merge in the astrophysically realistic context of cosmological 
structure form ation. In their analys is they take the analyti- 
cal formula of iLacev fc Coi3 l| 19931 ) and break it down into 
its principle components, fitting each parameter in order of 
importance using their measurements of mergers between 
0^z^2. J08 find that using alternate forms of Coulombs 
logarithm gives a better agreement in the mass ratio de- 
pendence. They also propose a new form of the circular- 
ity dependence that better accounts for radial orbits (where 
the formerly used p owerlaw breaks down). In contrast to 
iNavarro et al.l (Il995l). J08 find th at the commonly used for- 
mula from iLacev fc CoQ (| 19931 ) systematically underesti- 
mates the merger-time for minor mergers and overestimates 
the merger-time for major mergers. They conclude their 
study by providing a statistically robust and cosmologi- 
cally found e d form ula for merger-times. A follow up work 
IJiang et all l|2010t ) addresses minor issues in the baryonic 
physics such as overcooling. 

In the current work we extend J08s study in two ways. 
Firstly by confirming their result using a different area of 
structural evolution, i.e. looking at higher redshifts (z^lO) 
and a lower mass range. To this end we find our results are 
described well by J08s fitting formula. The second is by giv- 
ing an independent comparison with BK08 who published 



at around the same time. Both works look at the impact 
of orbital parameters on merger-timescales, however in very 
different ways. BK08 are able to make clean and accurate 
measurements as they are dealing with an isolated/closed 
system. J08 conversly, are looking at mergers amidst cosmo- 
logical structure formation where kinematics can be affected 
by the local enviroment. 

One difference between J08 and this study is the way 
in which merger-time is defined. We explore two types of 
merger; mergers that finish with coalescence of the two halo 
centres as well as mergers that are concluded by the tidal 
disruption of the satellite. These two merger criteria are ex- 
plained in greater detail in the method. 

It is necessary to highlight that the simulation used in 
this work contains dark matter only, while the study of J08 
used a hydro/N-body simulation including gas and star for- 
mation. Using a dark matter only simulation has limitations. 
With regard to mergers defined by halo cores coalescing, it is 
assumed that galaxies lie in the centre of their dark matter 
halos and while this is a reasonable assumption it does not 
represent a 1-1 mapping when studying effects such as AGN 
activity and SMBH growth. With respect to the tidal dis- 
ruption merger definition, the disruption of the dark matter 
halo does not necessarily correspond to tidal disruption of 
the galaxy. It is reasonable to assume that there is a connec- 
tion between the two and at the very least the disruption of 
the DMH serves as proxy for the tidal stripping time of the 
galaxy. 

The aim of this paper is to assess the validity of cur- 
rently used merger-time formulae with as many mergers (of 
sufficient resolution) as possible. The focus is to better ap- 
proximate the merger-times of DMHs in an astrophysically 
realistic context even if this does not correspond to an im- 
provement in the description of dynamical friction in ideal- 
ized isolated cases. 

We describe our method in section 2, present the results 
of our comparison in section 3 followed by conclusions in 
section 4. 



2 METHOD 

We ran a cosmological N-body simulation using the code 
GADGET2 iSpringcl 200^). We then identify DMHs using 
P-groupfinder. For each halo at z=0 we construct a merger- 
tree, which gives us a detailed dynamical history of all the 
halos in our simulation at all redshifts and allows us to look 
at all the mergers in the volume. Once mergers are identi- 
fied we track them in subsequent snapshots on a particle by 
particle basis. In the following sub-sections we describe the 
details of the above steps. 



2.1 Simulation and Group Finding 

A cosmological N-body simulation of dark matter was used 
to simulate a volume of the universe, 10 /i~^Mpc^ on one 
side, from redshift z=50 to z=0. The relevant cosmological 
parameters were fro m WMAP3 data: r2 M=0.24, f2A=0.76, 
0-8=0.74 and h=0.73 l|Spergel et al.ll2007l ). Our mass resolu- 
tion is 4.9x IO^Mq/i"^ (corresponding to 512^ particles) and 
we saved 102 snapshots between redshifts z=19 and z=0. 
The initial conditions (ICs) at z=50 were calculated 
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using the Zel'Dovich approximation, which assumes a hn- 
ear e volution in densi ty from z=1000 to the starting red- 
shift (IZerDovichlll97(]l ). It is true that using a higher start- 
ing redshift and more sophisticated methods of initiahzing 
the volume will res ult is a slightly different mass function 
(jCrocce et al.ll2006l ). we argue that this will not affect the 
measurement of merger-times for individual mergers. Anal- 
ysis of 512^ runs for 30 and 100 /i~^Mpc'^ volumes (with 
WMAP 5 and 7 year releases) is currently underway. The 
new volumes will be evolved from a redshift greater then 
250 and will allow a much more detailed description of the 
merging populations. 

We continue by identifying dark matter halos in a post 
simulation analysis with P-groupfinder using the Friends-Of- 
Friends (FOF) approach. We used the typical linking length 
of 0.2 times the mean inter-particle separation (correspond- 
ing to 3.9 kpc). 



2.2 Merger-Tree 

For every halo at z=0 we created a merger-tree. Our merger- 
trees were constructed in the conventional fashion whereby 
we take a parent halo from snapshot i and look for its largest 
contributor in snapshot i-1. Halos are linked to their ances- 
tors and descendants by unique particle ids. Some common 
pitfalls associated with building a merger-tree in this way 
are what we refer to as flybys and bridged halos. Flybys 
occur when the satellite halo enters the host halo only to 
continue on a perturbed (but ultimately unbound) trajec- 
tory in subsequent snapshots. Bridging happens when two 
small satellites appear to merge prior to falling into a much 
larger potential/halo. This is registered in the merger tree 
as two hierarchical mergers, but in reality the larger halo 
has simply accreted two small satellites. 

Both of these problems arise from the discrete nature 
of simulations outputs. Addressing the impact of these sys- 
tematic error s on galaxy evolution and black hole growth 
models (e.g. IVolonteri et al.l l2003l: ISoinerville et all l2008l: 



Micic et al.l2007l . l201ll : lKauffmann et al.lll993l : ICroton et all 
20061 ) is beyond the scope of this paper, however further re- 
search on this topic is important to understand and quantify 
the conseque nces this will have on current mode ls (Poole in 
prep, also see ISinha fc Hollev-Bockelmannll2012l '). In the ap- 
pendix we discuss some of the steps taken in the construction 
of our merger-trees to avoid false or artificial mergers. 



2.3 Halo Definitions and Merger- Tree Pruning 

We track halos with at least 500 particles. While the FOF 
halo mass is used to discern relevant halos in the merger- 
tree, the virial mass and radius are used in all other calcu- 
lations. We approximate the viral mass and radius by the 
spherically averaged over-density. We define a halo where 
pvir(z) = 200 pciit{z), in which pvir{z) is the average density 
within the virial radius and pcrit(z) is the critical density of 
the universe at the given redshift. Our merger-tree is con- 
structed using halo masses defined by the FOF-groupfinding 
algorithm. After applying the SO algorithm to these halos 
we find several of them lack a significant over-density (that 
would be classified as a core). These halos are naturally ex- 
cluded. 



We calculate the c entre of each halo potential using an 
iterative approach (e.g. lPower et al]|2003l ). We calculate the 
centre of mass of the group and proceed to remove particles 
beyond 0.98 Vi, where ri is the current radius (ro = r^ir). We 
then re-calculate the centre of mass for new set of particles 
and repeat the process until the number of particles reaches 
5% of the original set. Halo velocities are approximated using 
the weighted mean, where weights are based on the distance 
from the centre of the potential. 

We made several cuts to our collection of mergers in or- 
der to make our analysis set. As previously stated, there are 
a notable number of artificial mergers that we excluded from 
constructing our merger-tree. Flybys and bridged mergers 
have a significant impact on the accuracy of merger-trees, 
and we developed measures to exclude these from the analy- 
sis (s ee a p pend ix). In general phase-space merger-tree codes 
(e.g. iBeiiroozi et al.ll20ll] ) mitigate many of these artificial 
mergers automatically. Naturally we exclude any merger 
that does not finish by the end of our simulation. The effects 
of this can be seen as a distinctive curve in the upper left 
corner of the raw data plots of figure [T] This introduces a 
selection bias to our data set which will be discussed later. 



2.4 Defining Merger-Time 

Merger-time can be defined in multiple ways. Most involve 
the evolution of a characteristic or property of the merger 
such as separation, specific angular momentum or the num- 
ber of bound particles, each has its own pitfall. Separation 
becomes an inaccurate proxy for highly radial orbits. Tidal 
disruption can impede the measurement of kinematics when 
using a specific angular momentum criteria. The number of 
bound particles is a relatively robust condition but is best 
tuned to track mergers that entirely disrupt the satellite. 
It is, however, highly dependent on the criteria for disrup- 
tion and the force and mass resolution of the simulation 
iKlvpin et al.lll999h . 

It is important that we clarify the definition of a merger 
in its astrophysical context. In this work we use two def- 
initions of a merger, h owever we point the keen reader to 
IWetzel fc White! (2OI0I ) Fl. The types of merger assessed in 
this work are: (1) Coalescence, where halo cores (and their 
hosted galaxies) coalesce into a new potential before the 
satellite is tidally disrupted. This is t he classically pe r ceived 
pict ure of mergers in papers such as iTremaine et al.l (|l975t ) 
and IWhitd l|l976f ). (2) Disruption, where the satellite is 
tidally disrupted before the halo cores coalesce. The various 
stellar streams in our own galaxy (Bclokurov ct al. 2006 ) are 
testament to the ability of tidal forces to disrupt satellites 
before they coalesce with their host. 

The start of the merger is defined as the moment the 
satellite halo crosses the hosts virial radius. To get a bet- 
ter approximation of the exact moment we interpolate to 



1 IWetzel fc Wliit j 1I2OIOI ') define 4 types of satellite removal which 
are directed toward the abundance matching approach to galaxy 
evolution. Semi-analytical and hybrid models should intrinsically 
incorporate their 3rd and 4th types of removal if flybys and 
bridged mergers are dealt with in an appropriate fashion by the 
merger-tree. 
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first order all merger properties between the adjacent snap- 
shots. The end of the merger is quantified by two criteria: the 
evolution of specific angular momentum and the evolution 
of the number of bound particles. The Angular Momentum 
Criteria (AMC) is satisfied when specific angular momentum 
decreases to 5% of its initial value at the start of the merger. 
The angular momentum approach acts as a good criteria for 
coalescence and is the same approach taken by BK08. J08s 
halos have a baryons which can be used to directly track 
exactly when the two galaxies collide. The Bound Particle 
Criteria (BPC hereafter) covers tidal disruption. Similarly, 
the BPC is satisfied when the number of bound particles 
decreases to 5% of the initial value. We incorporate both 
types in a mixed definition, whereby the end of the merger 
is marked by whichever of the two criteria is satisfied first. 
When we compare our work to the literature we take care 
to note and match the merger criteria. 

2.5 Tracking Mergers 

We follow a fixed set of particles for the host and satellite 
halos at the time of in-fall. For the satellite we continue 
to use the same set of particles and retrieve their updated 
positions and velocities in subsequent snapshots. For the 
host, in order to take into account accretion, we use the 
original set of particles to track the central position and 
velocity of the halo but re-calculate the mass and radius 
based on the SO algorithm for all surrounding particles. 

This approach eliminates errors typically associated 
with tracking halos using group-finders. One such error oc- 
curs when two over-densities of comparable size occupy the 
same FOF halo. Then the group-finder tracks the largest 
of the two as the main halo, but when the over-densities 
are of comparable size the tag can artificially switch to th e 
other halo in the next snapshot (see lWetzel fc Whitell2O10l ). 
Another problem intrinsic to some subhalo finders (and any 
technique that relies on a density contrast to track substruc- 
ture) is that they can intermittently lose subs tructure when 
it pa sses by the dense centre of the host fsee iLudlow et al] 
I2OO9I ). The appendix highlights the benefit of using a fixed 
particle set to track mergers. 



3 RESULTS 

Throughout the results section we compare the difi'erent 
types of end-of-merger criteria and how they vary with or- 
bital parameters. Recall that in both the AMC and BPC 
criteria the merger is considered complete when the value is 
5% of its original value at in-fall. When we use AMC only, 
our data set consists of 1116 mergers (denoted by red un- 
less stated otherwise). When BPC is applied the data set 
consists of 1199 mergers (blue). In the mixed set the end 
of the merger is defined by whichever criterion is satisfied 
first. This gives us 1438 mergers (black). Interestingly, 877 
of the mergers have end-time measurements for both AMC 
and BPC methods. Furthermore, 55% of this sub-set have 
Tbpc < Tamc- Throughout the results section any data 
bin used to compare variables contains a minimum of 40 
mergers. The only bins that do not have N ^ 40 are distri- 
bution functions, usually in the bottom panels. Overall we 
find that merger-time is most strongly correlated with the 



dynamical timescale of the host halo. With regard to other 
merger parameters we get varied results for the different 
merger definitions. 

The results are presented in the following form; in sec- 
tion 3.1 we look at the merger-time vs dynamical time re- 
lation, in section 3.2 we investigate mass ratio, circularity 
and orbital energy and finally in section 3.3 we compare our 
measured merger-times with the predictions of other merger- 
time formulae and provide our own best fit parameters. 

3.1 Merger-Time with Redshift 

Merger-times for such high redshift halos has remained un- 
explored in cosmological simulations. Figure [T] shows the 
measured merger-time as a function of the redshift at which 
the merger starts. The left panel shows data from the AMC 
and BPC merger sets separately. The central panel shows 
the mixed data set while the right panel shows the binned 
equivalent of each merger set. Bins are spaced equally in 
starting redshift. Points represent the bin median and error 
bars show the inter-quartile range. The red curve shows the 
limit of our measurements imposed by the end of the simu- 
lation at redshift z=0 (if the merger does not end in Hubble 
time, it is excluded from the data set). The measurement 
limits introduced by the end of the simulation will give rise 
to a selection bias for our low redshift mergers. This can be 
seen in the first couple of bins in the right panel in figure[T]as 
data starts to deviate from the fit. In general, all three sets 
of data show that high redshift mergers take considerably 
less time to finish then those at low redshift 

One of the principal components of most prescriptions 
of merger-time is the dynamical timescale of the host halo 
Tdyn. In the context of halo mergers, dynamical time is the 
time necessary for a satellite halo to make its first pericentric 
passage of the host halo. The redshift evolution of a halos 
dynamical timescale in spherical collapse models can take 
multiple forms, 

T,„„ = -I^ = o.m(^)-\ (2) 

where H{z) is Hubble's parameter and for the relevant cos- 
mology is given by H^{z) — Hq{Qa + O.m.a'^)- Note we leave 
out a factor of -\/2 in equation 2 as this will later be absorbed 
in other constants. The growth of a halos mass with redshift 
enlarges the virial radius as well as increases the dynami- 
cal time. Both BK08 and J08 take this effect into account. 
BK08 fit their orbital parameters to Tmergo/Tdyn (where Tdyn 
is constant for their host) while J08 have a unique r^ir/Vc 
for each merger. The right panel in figure [1] shows how the 
merger-times measured in our simulation correlate with dy- 
namical time. We plot 3 multiples of Tmerge = aTdyn to 
illustrate the impact of Tdyn in our results. We do see indi- 
cations that there is a redshift dependence beyond the Tdyn 
correlation. Such a dependence might manifest in the red- 
shift evolutio n of the distribution of orbital parameters (like 
IWetzelll201ll ). We hope to study this in greater detail with 
the new simulations. 

3.2 Other Orbital Parameters 

The dependence of DMH merger-times is often tied to other 
orbital parameters such as mass ratio, circularity of the orbit 
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Figure 1. Merger-time against starting redshift. The plot holds a total 1438 mergers with a mixed end definition (threshold 0.05). The 
left panel shows the measurements for the AMC and BPC merger sets along with the associated errors. The middle panel shows the 
data for the mixed merger definition. The red curve for the first two panels represents the maximum merger-time measurable at a given 
starting redshift. The plot on the right shows the data for all sets binned in starting redshift. The central point for each bin corresponds 
to the median while error bars represent the interquartile range. 



20 



15 



10 



Fit /3 = 1-3 7 = 1.0 

Fit (3 = 1.6 7 = 0.7 

Fit /3 = 2.4 7 = 0.6 

AMC ( 1 I I 5 mergers) 

Mixed (1438 mergers) 

BPC (1 199 mergers) 




Figure 2. A comparison of Tincrgo/Tdyn vs mass ratio for the 
different end criteria. The central point in each bin represents the 
median while error bars correspond to the interquartile range. 
Each set is fitted to the equation shown. The bottom panel shows 
the fractional distribution of the merger sets with mass ratio. The 
fit values are shown in the top right. 



and the orbital energy. In this section we explore these pa- 
rameters and their dependence on the length of the merger. 
Theory, especially with respect to the classically per- 



ceived picture of mergers, expects a rather strong depen- 
dence of merger-time on mass ratio. There have been a num- 
ber of works aiming to refine the form of the Coulomb log- 
arithm. It is becoming increasingly accepted that the mass 
ratio dependence of merger-time is best characterized by: 

Mp/Ms 



T, 



ln(l + Mp/M,)' 



(3) 



where the Coulomb logarithm is ln(l + Mp/Ms). Figure [2] 
shows measured merger-time (corrected for the Tdyn depen- 
dence we see in figure [1} as a function of the mass ra- 
tio of merging halos. In all sets of data there is a clear 
trend that the duration of the merger increases as the 
host halo becomes larger than the satellite. Our results 
for the AMC mergers are consistent with previous works 
jWetzel fc WhitelboiOl ). While there is considerable scatter, 
the median value for each bin is in good agreement with 
the fit. The general scatter can be attributed to freedom in 
other orbital parameters. Despite fitting to the functional 
from BK08 (shown in figure [2} we find an exact match with 
JOS, i.e that 7 = 1.0. 

When the BPC is used, we see that merger-time does 
not vary with mass ratio as expected. While the trend is 
correct, the dependence is weaker than seen in other works. 
The histogram in the bottom panel of figure [2] shows the 
fractional distribution of merger sets with respect to mass 
ratio. 

In any closed two body system the exact trajectory can 
be described by circularitjo and orbital energy. Both of these 
variables are hard to constrain in such a broad parameter 
space. As seen in JOS it is sometimes easier to look at weaker 
dependences by plotting them against Tsim/Tmodci. 



Circularity (e) is related to the eccentricity(e) by e = \/l — e^. 
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Figure 3. Circularity as a function of orbital energy. The top 
panel shows raw data with the selection bias (see text) marked by 
the red line. Also shown in the top panel is the unbiased samples 
used for the frequency distributions in the lower two panels. 



Figure[3]shows orbital energy as a function of circularity 
for all mergers used in our analysis (i.e the mixed set). A 
selection bias can be seen in the top panel in the form of 
an 'envelope' of points. This is the product of the definition 
of merger start time (i.e the point at which the satellite 
crosses the hosts rvir). The detection of mergers with higher 
orbital energies and/or more circular orbits is limited by 
'"per/f'vir !^ 1; in Other words we will only detect mergers 
with orbital pericentres inside the hosts virial radius. An 
orbits pericentre is defined e). For circular 

orbits the semi- major axis a is equal to rc- This constraint 
can be expressed in terms of the circularity of the orbit such 
that 
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n 
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F(e) = 1.4exp[0.5E] Mixed (^438 mergers) )K 
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Figure 4. Merger-time as a function of circularity. Points cor- 
respond to the median while errors represent the interquartile 
range. The bottom panel shows the fractional distribution of the 



merger sets with respect to circularity, 
bound and mixed particle sets. 



X 



fits are shown for the 



rc(l-e) 



1. 



(4) 



In the limit of this constraint rpcr/7"vir = 1, equation U can 
be re-arranged to give, 



1 + 



(5) 



Equation [5] is shown in figure |3] as a red curve. It is in ex- 
cellent agreement with the envelope. 

The bottom two panels of figure [3] show the distribu- 
tion of circularity and orbital energy respectively. In each 
panel we show the distribution of the full data set as well as 
an unbiased sample. Sub-samples are set by equation [S] and 
highlighted with lined boundaries in the top panel. The mid- 
dle panel shows that selection bias has a strong influence on 
the distribution of circularities. This bias will not influence 
the fit of merger-time (section 3.3) but may have an impact 
on the fitting of circula rity in f i gure Th e unbias e d set is 
in good agreement with IWetzel (2OIII ') and iBensoiil (120051 ). 

The dependence of merger-time on circularity is less 
constrained than the mass ratio relation. Figure U shows 
merger-time against circularity for all 3 merger sets, cor- 
rected for the already fitted parameters. The overall trend 
shows that the duration of mergers increases with increasing 
circularity. 

This positive correlation of merger-time and circularity 
is seen most strongly in the BPC set. Under the BPC defini- 
tion, the trend corresponds to radial orbits tidally disrupting 
faster than circular orbits. This is expected as radial orbits 
make close central passings where tidal forces are strongest. 
The AMC merger set deviates from the trend significantly 
for more radial orbits. This deviation persists when an unbi- 
ased sample of circularities is usecjfl- It appears to be physical 



^ It also persists under the condition that the satellite has a well 
resolved core. 
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and not an artifact of method or numerics. The imphcation 
for merger models is that, for well resolved satellites, both 
highly radial and highly circular orbits take a similar time 
to coalesce. Arguably, the binned points in figure 9 of JOS 
show a similar trend. Their trend is much weaker and it 
would require a detailed analysis to confirm it is a physical 
effect. Radial orbits e < 0.3 are often overlooked in studies 
of this type, so perhaps this effect has not been appreciated. 

We fit the mixed and BPC merger sets to the depen- 
dences seen in JOS and BKOS (shown top left ). JOS have re- 
tained the analytically founded forirQ seen in lLacev fc Cold 
||1993D while BKOS use an imperically motived exponent. 
While both forms of the dependence are in accord with the 
data, the model with an extra fitting parameter (JOS) shows 
a tighter agreement. We are not able to fit the AMC merger 
set due to the upward trend at low circularities. The bot- 
tom panel of figure |3] shows the frequency distribution of 
each set. 

The final theoretical dependence is the orbital energy of 
the merger. As previously stated the proxy for the energy of 
the orbit is rc{E) /rvir- This compares the radius of a circular 
orbit (with the same energy as the orbit in question) with the 
virial radius of the host. The top panel of figure [S] compares 
the dependence of rc{E)/rvii with merger-time. Similar to 
previous figures, other dependencies on merger-time have 
been taken out. The exception is the circularity dependence 
for the AMC set. Overall there is a weak dependence in 
all 3 data sets. Other works suggest Tmerge oc [rc(i?)/rvir]'', 
where /x ^ 1. BKOS find fj, = 1. JOS do not include an orbital 
energy dependence in their equation 5 for simplicity and ease 
of use. They do however propose an alternate form of the 
equation in the conclusion that incorperates the satellites 
orbital energy. We find a relatively consistent value of /i = 
0.1 across all data sets. 

The fractional distribution of rc{E)/rvii is displayed 
in a histogram in the lower panel. The range of rc(iJ)/rvir 
found in our simulation extends beyond the typically shown 
range. The peak of the distribution is slightly beyond 1. A 
close examination (snapshot to snapshot) of these extremely 
energetic orbits usually shows a close passing of a third body 
also orbiting hos t halo (a relatively common occurrence e.g. 
ISales et al.|[2007l ). 

3.3 Comparison with previous works 

In this section wc compare our measure ments with the em - 
piri cally derived formulae of |Bpvlan- K olchin et al.l (|200Sl ) 
and I Jiang et~all (20081 ) . We remind the reader that the com- 
parison is only with measurements from the AMC merger 
set as this is the merger definition under which their formu- 
lae were constructed. The best way to make this comparison 
is to look at the variance of the ratio of measured/predicted 
merger-time (Tsim/Tmodci) for different orbital properties. 
The formulas used in our comparison are 



\Bovlan-Kolchin et al\ \200^) : 



■ tdy 



AjMh/Ms)'' 
Tn(l + Mh/M, 



-exp 



C 



ME) 



rc{E) 



(6) 



3.0 



2.5 
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0.5 

0.0 
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Fit /c=1.19, /i = 0.1, (not corrected for F[£]) 
Fit ic = 0.87, p. = OA 
Fit /c = 0.90, u = 0,1 



F[r,(E)/rJ = /c(r,(E)/r J" 
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AMC (1116 mergers) □ 
Mixed (1438 mergers) ))( 
BPC (1 199 mergers) A 




rc(E)/r„ir 

Figure 5. The top panel shows merger-time against rc{E) /r^i^ (a 
proxy for orbital energy). The bins are spaced equally in log. The 
distribution is skewed towards lower ratios, as shown in the lower 
panel. The fitted parameters are shown top left. Points corre- 
spond to the median while errors show the interquartile range. 



where J/Jc{E) = e and their fitted constants are 
A = 0.216, B = 1.3, C = 1.9, D = 1.0. 

\ Jiang et al] \200^) : 



Trr 



f{e)Mh rvir 
2CMs InA Vc{vir) 



(7) 



(8) 



where /(e)and InA are the corresponding fitted values from 
their worC|, i.e /(e) = 0.94e''-^ -I- 0.60 and our eqE] 

In addition to the above formulae, we do a maximum 
likelihood estimate of equation[S]using our merger-time mea- 
surements and find the following values. 



A = 0.9, B = 1.0, C = 0.6, D = 0.1. 



(9) 



We use these parameters to compare our fit with the other 
two works in figure [S] It shows the variance in Taim/T^odei 
with mass ratio, circularity and orbital energy. From top to 
bottom, figure |6] shows a comparison of BKOS, JOS and the 
fitted values in this work. BKOS formalism deviates from our 
measurements as a function of mass ratio and orbital energy. 
At high host-to-satellite mass ratios BKOS overestimate the 
merger-time by a factor of ~ 3. The power law in the mass 
ratio dependence of BKOS is responsible for the deviation. 



* Introducing an additional fitting parameter to accomodate the 
power law for highly radial orbits. 



^ We also compare with the alternate form of their fit (their equa- 
tion 8) but find it to be a worse fit. 
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Figure 6. Tsim/Tmodel against mass ratio, circularity and orbital energy. In this plot we compare the empirically derived formulae of 
BK08 and JOS to the parameters found in this work. Points correspond to the median while errors show the interquartile range. 



The discrepancy between BK08 and our fitted parame- 
ters (eq [7] and (Qjl indicates that the dependence of merger- 
time on circularity and orbital energy is weaker than pre- 
viously derived. This illustrates the difference between an 
idealized isolated merger simulation and a full cosmological 
simulation. BKOSs mergers were in isolation and had relaxed 
Hernquist halos. Conversely, the mergers in this work are in 
a cosmological context, here the host is still collapsing, un- 
dergoing multiple mergers at once and at the mercy of the 
tidal fields in its local environment. 

BK08 explicitly state the range of validity of their 
formula as 0.025 Ms/Mp < 0.3, 0.3 < e ^ 1 and 
0.65 rc{E)/rviT ^ 1 . We find that the range of valid or- 
bital parameters corresponds to a significantly limited range 
of mergers in our simulation. A comparison within this pa- 
rameter space can be seen in figure[71 It shows a very modest 
number (~ 3% of our AMC set) of comparable mergers. 

The comparison to JOS produces a much tighter agree- 
ment in all three parameters (figure|6] middle panels). Max- 
imum deviation occurs for equal mass ratio mergers as well 
as highly radial and energetic orbits. Our own fit to BKOSs 
functional form shows a slight reduction in scatter when 



compared with JOS. This may just correspond to the differ- 
ent ways the mergers are qualified. Explicitly that we use 
specific angular momentum while JOS are able to measure 
when the baryons in the galaxies collide. Either way the 
results are very well described by the formula of JOS. 



As a note, caution should be taken when applying 
merger-time formulae. The accuracy of a formula can be 
highly dependent on the conditions under which it was de- 
rived. For example, as the merger progresses from some fidu- 
cial starting point the parameters of the merger change. The 
mass of the satellite typically decreases while the mass of 
the host increases. rc{E)/rvii tends toward zero, due to the 
growth of the virial radius with cosmic expansion as well 
as the dissipation of orbital energy from particle-particle 
interactions. These inputs will also vary depending on the 
halo finder/halo definition used to construct the merger-tree. 
The application of any approximating merger-time formal- 
ism should be consistent with the circumstances under which 
it was derived. 
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Figure 7. A comparison with BK08 showing the predicted and 
measured merger-times for mergers within their valid parameter 
range. 



4 CONCLUSIONS 

This paper presents results of a cosmological N-body sim- 
ulation in a small cosmological volume where dominantly 
spiral and dwarf galaxies form. We studied the mergers of 
the dark matter counterparts at high redshift. The mass and 
redshift range of our simulation has never been probed for 
merger-times. Our study considers two merger definitions 
that qualify the end of the merger, tidal disruption and core 
coalescence. We compare the measured merger-times from 
our simulation with the predicted times from other works. 

The merger-times of DMHs is strongly dependent on 
the redshift at which the merger starts. This manifests it- 
self in the redshift evolution of the host halos dynamical 
timescale. The dynamical timescale of a halo corresponds 
(approximately) to the time it takes for a point mass to 
make its first pericentric passage. 

We consider both the AMC and BPC merger definitions 
(see section 2.5) in our analysis. Both scenarios impact on 
satellite removal in N-body /semi-analytic recipes of galaxy 
evolution. From the mergers that have both AMC and BPC 
end measurements, we find that 55% have Tbpc < Tamc- 
We also find a selection bias that arises from the pre- 
requisite that mergers included in our analysis must have 
rpcr/rvir ^ 1- This bias is characterized by equation [5] and 
should be considered in kinematic works of this type. 

In the classical merger definition (AMC - where resolved 
satellites coalesce) we find good agreement with Tdyn and 
mass ratio relations. The circularity and orbital energy have 
less of an impact on merger-time. Our results show a much 
weaker dependence on orbital energy than analytical/ideal 
works. We find an interesting result with regard to the cir- 
cularity dependence; for well resolved halos, highly radial 
and highly circular orbits finish in comparable time. This 



low circularity region (e < 0.3) is seldom probed, so it is 
possible (although unlikely) that this effect has been over- 
looked. BK08 do not cover e < 0.3 in their study and JOS do 
not show as strong an upward trend for low circularities. We 
have taken steps to check for numerical and method based 
artifacts and plan to perform a more in-depth analysis in 
the future. For e > 0.3 we find reasonable agreement with 
other works. 

Our results according to the tidal-disruption merger- 
criteria (BPC) show some interesting variations. We found 
a tight correlation with the dynamical timescale of the host 
halo. We also see a dependence on mass ratio, although it is 
weaker then expected (see figure[3]). Circularity for the BPC 
set shows shorter mergers have more radial orbits. In such 
orbits the satellite is exposed to stronger tidal forces. 

We compare o u r me asured merger-tim es wit h 
iBovlan-Kolchin et al.l (|2008l ) and Ijiang et al.l (|2008h . 
We find a systematic deviation between our measurements 
and the predicted times of BK08 who numerically fit their 
formula from isolated mergers. By fitting our cosmological 
mergers to BKOSs formula (eq ^ we are able to directly 
compare idealized and cosmological merger-times. As 
such, we find the biggest difference is the dependence on 
orbital a power of 0.1 instead of 1). 

Our results also show a mass ratio power of 1.0 (formerly 
1.3) and a milder circularity dependence. We cannot say 
exactly why there is such a contrast between mergers in a 
cosomological context and those in isolation, however part 
of this discrepancy can be attributed to the fact that in 
cosmological structure formation multiple mergers occur 
simultaneously. In this work we have empirically found the 
differences between the two cases and the next step is to 
understand the physical mechanisms responsible. 

In contrast to BK08, our data is well described by J08s 
approximating formula. Their empirically fitted formula was 
derived under the same cosmological N-body environment 
as this work. The differences between the simulations in JOS 
and this work are the inclusion of baryons as well as the mass 
and redshift range covered. We remind the reader that there 
are limitations to approximating the AMC and BPC merger 
measurements using dark matter alone. A comparison of our 
fit with JOS shows a slight reduction in scatter, however this 
is consistent with the different approaches used to measure 
merger-time. While both formulas fit the data well, the for- 
mula of JOS has one less input parameter which makes it 
easier to use. Our fit to BKOSs form provides slightly re- 
duced scatter at the expense of a additional parameter. 

A comparison of the two works (isolated to cosmo- 
logical) illustrates that while isolated/idealized numerical 
experiments are important to understand the underlying 
physics, a more versatile (statistically robust) approxima- 
tion may serve better in non-ideal applications, i.e. semi- 
analytic/hybrid models. 
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6 APPENDIX 

6.1 Details of our Merger-tree 

ISinha fc HoUev-Bockelmannl (l2012i ) try to quantify the effect 
of flybys for corrections in current models. To address this 
issue and avoid "false mergers" in the current work we apply 
the following conditions in the construction of our merger- 
trees. 

• Any halo associated with the tree has at least 50% of 
its particles in the "root" halo at z=0. 

• If there is any ambiguity as to which halo is the main 
progenitor (in snap i-1), we then follow the contributions of 
the both progenitors back to i-10. 

• One of the benefits of P-groupfinder is that parti- 
cle ids in a given halo are listed in order of binding en- 
ergy (to a large extent although not exactljjf]) this allow 
us to weight the mutual particles such that more bound 
p artic les hold greater weight when considering progenitors 
(|Bovlan-Kolchin et al. 1120091 ) . 

• When a halo in snapshot i traces to a significantly larger 
halo in snapshot i-1 (i.e an emerging flyby), we look at snap- 
shots i-10 through to i-f 10 to firstly establish what the halos 
mass was upon in-fall and secondly to see whether the halo 
grows significantly after it emerges. To distinguish between 
flybys and highly eccentric orbits we employ the condition 
that a satellite must at least double its initia l in-fall mass to 
be considered independent after emergence (|Wetzelll201ll ). 

While these approaches help to minimize the number 
of flybys it does not catch all of them. A small fraction of 
mergers are false mergers due to flybys. These mergers are 
removed from further kinematic analysis. Figure [8] shows an 
example merger-tree. 

6.2 Our Tracking Scheme: The Importance of a 
Fixed Particle Set 

An initial attempt to track mergers by comparing substruc- 
ture at different snapshots (via common particle ids) proved 
difficult due to the limitations of single epoch group-finders. 

^ P-groupfindcr lists the particles in order of binding energy for 
each substructure. It lists the substructure in a given halo from 
largest to smallest. 
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Figure 8. An example of one of our merger-trees. The size of the 
circles is proportional to the size of the halo. 



Subfind identifies substructure by constructing a 
smoothed 3D density contour of the FOF halo. As the global 
density is lowe red, saddle points e merge between adjacent 
over-densities (|Springel et al.ll200lf ). This gives a subset of 
particles that Poisson's equation can be applied to. 

Subfind can temporarily loose substructure when the 
subhalo is passing through an over-dense region. The cause 
of this intermittent loss arises from the position based (i.e. 
smoothed density) selection of particles. It is important to 
track all particles associated with the satellite halo at in-fall 
when measuring the evolution of satellite properties. 

Figure shows a 2D projection of structure found by 
various methods at three points during the merger. The first 
two columns show the subfind substructure(blue) compared 
with the satellite tracing scheme used in this work (green). It 
is important to highlight that the blue particles correspond 
to all substructure in the FOF halo and the plots in a each 
row are centred on the same coordinates. The middle row 
(z=2.22) shows the substructure found as the satellite passes 
through a region of high background density. 



The third column of figure [9] has a different scale to 
the first two in order to highlight the bound particles in 
reference to their unbound counterparts. This demonstrates 
the effectiveness of our tracking scheme above a nd beyond 
the gr oupfinder-particle tracking approach. See iHan et al] 
for a more detailed comparison. 
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XY plots of a satellite halo during merger 
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Figure 9. A 2D comparison of tiie structure found witii subfind at 3 points during a merger. Tlie first two columns sliow tiie subfind 
substructure (blue) compared with the satellite tracing scheme used in this work (green). The final column shows the bound(grccn) to 
unbound (red) particles of the satellite in the broader context of tidal stripping. Note that column 3 has a different and varying 
distance scale to the first two columns. All columns are centred on the same coordinates. 
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